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Abstract 

Accurate and sensitive measurements of higher order cu- 
mulants open up new approaches to pion interferometry. It 
is now possible to test whether a given theoretical prediction 
can consistently match cumulants of both second and third or- 
der. Our consistency test utilizes a new technique combining 
theoretically predicted functions with experimentally deter- 
mined weights in a quasi-Monte Carlo approach. Testing a 
general quantum statistics-based framework of Bose-Einstein 
correlations with this technique, we find that predictions for 
third order cumulants differ significantly from UA1 data. This 
discrepancy may point the way to more detailed dynamical 
information. 

PACS: 13.85.Hd, 25.75.Gz, 12.40. Ee 



Pion interferometry has been a vital part of multipar- 
ticle physics for several decades (!]]. While traditionally 
experimental effort in this field has centered around sec- 
ond order correlations, much progress has been made 
recently in accurately quantifying higher order correla- 
tions by means of so-called correlation integrals pj, to 
the point where these now yield statistically significant 
conclusions not only for moments but also for higher or- 
der cumulants. Because cumulants are so sensitive to de- 
tails of the dynamics, they represent a stringent testing 
ground for proposed theoretical models. 

A number of theoretical predictions for higher orders 
exist p-pf. In particular, Andreev, Pliimer and Weiner 
(APW) |(J have suggested a very general quantum- 
statistical framework, based on the classical source cur- 
rent formalism applied successfully in quantum optics. 
Its basic assumptions are (1) a Gaussian density func- 
tional for the classical random currents and (2) isotropy 
in isospin space. These two assumptions determine all 
higher order correlation functions in terms of the basic 
correlator, independent of the structure of the sources. 
All further assumptions concern only a more detailed 



specification of the space-time evolution of the sources. 
Thus, the APW framework includes as special cases more 
specific models of Bose-Einstein correlations such as the 
GKW model and the approach of Biyajima et al. [Q. 
Because the APW model is so important, we test a sim- 
ple version of its predictions below. It will also serve as 
an example to show how our approach works. 

While higher order cumulant measurements are valu- 
able in their own right, they can be used to even greater 
effect in consistency checks: once an assumed parameter- 
ization is found to fit the second order data, the same set 
of parameter values ought to fit all predicted higher or- 
der correlations as well. Departing from tried and tested 
ways, we therefore concentrate not so much on numerical 
values of source parameters, but rather on utilizing their 
required constancy over cumulants of different orders to 
test for consistency and ultimately falsifiability of a given 
theoretical prediction. 

Pion interferometry measures correlations in terms of 
pair variables such as 3- or 4-momentum differences of 
two particles. In the correlation integral formalism M, 
this means that the original variables in the rth facto- 
rial moment density p r {pi, ■ ■ ■ ,p r ) must be converted to 
r(r — l)/2 (not necessarily independent) pair variables 
(912, 913, ■ ■ ■ , 5r-i,r)- Here, we use the specific definition 
of a "distance" \p t —pj\ = y/{pi — Pj) 2 - (E~i - Ej) 2 be- 
tween two phase-space points pi and pj, but other pair 
variables and distance definitions are also possible. The 
differential moment densities as well as the correspond- 
ing cumulant densities can be written in terms of delta 
functions involving pair variables, byi = S(qi2 ~\pi—P2\)- 
For example, in third order we have 



P3(9l2,g23,93l) = Jdpidp 2 dp 3 P3(pi,P2,P3)8l26 2 3d; 



31 



The pair-variable moments are measured by direct and 
exact sampling prescriptions 0] in terms of the measured 
four-momentum differences between track pairs, Qij, and 
event averages (• • •), 



P2 = (Ei^j %12 -Qij)), 



(1) 



P3 = (Ei^fc %12 - Q»i)%23 - Qjk)S{q 3 i - Qki)) , (2) 
etc. Pair-variable cumulants 

£2(912) = P2(Ql2) ~ Pi®Pi(Qi2) , (3) 

C3 (912, 923, 93l) = PS ~ L(3)/°2®Pl + 2pi<gipi®pi , (4) 

are similarly measured in exact prescriptions. Mixed 
terms such as p\®p\ require multiple event averages in 
terms of distances Qfj = \pf — p b , | between particles pf 



and p b taken from different events a and b, 



(5) 



and so on (see Refs. g for details). 

These experimentally measured cumulants are, after 
normalization, to be compared to the corresponding the- 
oretical predictions. The APW normalized cumulant pre- 
dictions are built up from normalized correlators dij, 
the on-shell Fourier transforms of the space-time clas- 
sical current correlators. The specific parameterizations 
of the correlator we shall be testing are, in terms of the 
4- momentum difference qa, 



Gaussian: dij = cxp(— r 9y) , 
exponential: dij = exp(— rqij) , 
power law: dij = q^ a . 



(6) 

(7) 
(8) 



For constant chaoticity A and real-valued currents, APW 
predict the second and third-order normalized cumulants 
&2 h (9i2) and k^(q 12 , 923, 931 ) to have the form [|| 



th _ 
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piC 



= 2A(1 - A)di2 + \ 2 d\ 2 , (9) 

= 2A 2 (1 - A)[di 2 (i23 + ^23^31 + d 31 d 12 ] 

+2A 3 d 12 d2 3 d 3 i , (10) 



with any of the above dij parameterizations. 

A comparison of sample cumulants (0)-(H) with the 
theoretical APW cumulants in (p|)-(|lO|) encounters two 
basic difficulties: projection and normalization. In deal- 
ing with these, we are led to the "Monte Karli" prescrip- 
tion outlined below. 

Statistical sampling limitations render a fully mul- 
tidimensional measurement of cumulants such as 
£3(912,923,931) impossible, necessitating projections 
onto a single variable. Such projections must of course 
then be applied to theoretical predictions such as (|10| ) 
also. However, the complicated combinatorics of dij 's en- 
tering the cumulants make it generally impossible to give 
analytical formulae in terms of any one combined vari- 
able (such as S = <7i2 +923 + 931) without making further 
approximations. While two such approximations, leading 
to simple third order formulae, have been used previously 
by us and other groups |q-|lC|] , conclusions based on them 



are weak because the effects of making such approxima- 
tions are seldom quantifiable. 

The difficulty in normalization arises because ex- 
perimental and theoretical procedures differ. Theo- 
retical cumulants are normalized fully differentially 
with pi<g>Pi(9i2) and pi<8>pi®pi(9i2,923,93i) respectively. 
Such normalizations assume essentially perfect measure- 
ment accuracy for the qij 's and infinite sample size. Ex- 
perimentally, however, one always measures unnormal- 
ized cumulants over some phase-space region VL and then 
divides by an uncorrelated reference sample integrated 
over the same $7. For example, in second order one actu- 
ally measures 



A AT (^ 



kCni)dq 
J n Pi®pT(q)dq ' 



(11) 



which approaches the differential k\ {q) or the bin- 
integrated normalized theoretical cumulant 



Atff(fi) 



« 2 



(9) dq 



(12) 



only in the limit fi — ► and may significantly differ from 
AKf 1 for larger fi. Non-stationary single-particle distri- 
butions can significantly distort results, especially when 
they vary strongly within a given bin, as they do in the 
case of momentum differences. The problem is exacer- 
bated in higher orders. 

"Monte Karli" (MK) integration provides a solution 
to both these problems |g|. It is conceptually quite sim- 
ple: Since 6*2(9) = ^2(9) Pi®Pi(q) we can, in analogy to 
AA£ X , define an MK normalized cumulant 



M MK (fl) s Inkf{q)p l ®pT{q)dq 



f n pi<»pf c (q)dq 



kCfH<l)dq 
f n p-L<8>pl*(q)dq ' 

(13) 



where kf 1 is a function supplied by theory while pi®p1 x 
is taken directly from experiment. Thus AK 2 V[K by con- 
struction satisfies both the theoretical function k^ 1 and 
the experimental one-particle distributions. Drawing k 2 l 
into the event averages of the correlation integral pre- 
scription (0), the unnormalized MK cumulant becomes 



c^(q) = ((E l .A<i-Qt)K h (Qt))), 



(14) 



so that the argument of the theoretical k^ is now an ex- 
perimental particle pair distance, obtained by event mix- 
ing. Prescriptions such as (111) and dl^) below amount 
to Monte Carlo-style sampling of the theoretical func- 
tion weighted by the actual experimental particle distri- 
butions, thereby earning Monte Karli its name. Inte- 
grated over the bin domain f2, the normalized MK cu- 
mulant reads 



AAT K (tt) 



({J2i,jJadlS(q-Qtm 



(15) 



where binning of track pairs is represented by J„ dqS(q — 
Qfj): whenever a pair's Qf & falls within il, the numerator 
counts k^iQij) while the denominator counts 1. 

Similarly, the unnormalized third-order Monte Karli 
cumulant C^ K (q) can be written in terms of kf 1 times 
P\®pi®Pi^- For the projection onto the single variable 
q, different "topologies" exist which quantify the size of a 
given cluster of r particles in phase space Q . The "GHP 
max" prescription used here bins cluster sizes according 
to the maximum 4-momentum difference. Inserting this 
into the correlation integral prescription, one finds 



C? K (q) = ((j25[q-™*(Q 



k th (n ab n bc t 



ij ' ^°cjki ' 
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(16) 



logarithmically between q = 1 GeV and 30 MeVQ. Fits 
to the data were performed using the three parameter- 
izations (§)-(|) in the APW form (|). All fits shown 
include, besides the free parameters A and r (or a), an 
overall additive constant A as free parameter. Such addi- 
tive constants are necessary for our choice of normaliza- 
tion because of the non-Poissonian nature of UAl data, 
which leads to nonvanishing cumulants at large q. Best- 
fit parameter values obtained were A = 0.050±0.005, a — 
0.64 ± 0.03 (A = 0.269 ± 0.012) for the APW power law, 
and A = 0.59±0.04, r = 1.16±0.02 fm (A = 0.385±0.003) 
for the APW exponential. The Gaussian parameters are 
A = 0.23 ± 0.01, r = 0.81 ± 0.02 fm (A = 0.396 ± 0.002). 
Goodness-of-fits were x 2 / Nr > F = 3.34, 6.65 and 20.0 for 
APW power, exponential and Gaussian respectively. 



and the normalized binned MK cumulant becomes 



(((£i,j,fc/fi^<%-max(- 



■)]kf{Qll 
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(17) 
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Multiple event averages are handled efficiently within a 
"reduced event mixing" algorithm using unbiased esti- 
mators, cf. Eggers and Lipa in |g]. Because the identical 
projection is applied to both experimental and theoret- 
ical cumulants, their comparison is not subject to any 
approximation or limitation. Clearly, the MK procedure 
can be generalized to all possible moments and cumu- 
lants, independently of variable or integration topology. 

We have measured second- and third-order normalized 
cumulants using a sample of around 160,000 non-single- 
diffractive events at 630 GeV/c recorded in the UAl cen- 
tral detector. Only vertex-associated charged tracks with 
transverse momentum p± > 0.15 GeV/c, pseudorapid- 
ity 1 77 1 < 3 and good measurement quality were used. 
We restricted our analysis to the azimuthal angle region 
45° < \<f>\ < 135° because there acceptance corrections 
are small enough to be safely neglected. To calculate the 
four momentum difference q, all particles were assumed 
to be pions. The resolution Aq ~ 8 MeV, estimated 
from the errors of track fits, is approximately constant 
over the whole range q < 1 GeV. For further details of 
the detector and pair selection criteria, see Ref. | llf| . 

In the following, we check the APW formulae ( 6^— (|10|) 
for consistency with the data. The procedure is to find 
numerical values for r (or a) and A by fitting the model's 
second-order prediction (0) for fe to the experimentally 
measured differential cumulant. These best-fit values are 
ported into the third-order prediction (J1C| ) which are then 
compared to third-order data using (|17|). 

In Figure 1, we show the second-order cumulant of 
like-sign particles AK2 = {J P2/ j Pi®pi) — 1, where nu- 
merator and denominator are integrals over bins spaced 



(++) and (-) 



j UAl data 

— APW power law 
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FIG. 1. Second order like-sign UAl cumulant with fits us- 
ing various parameterizations for di2 in formula (ffl). 

For AK3, we have performed three separate consis- 
tency checks, two based on approximations (shown in 
H ) , the third involving the Monte Karli method (shown 
below) which is exact. In Ref. S, the APW predictions 
for AK3 based on the two approximations were found to 
differ substantially from the data. 



1 It should be remarked that previous work has shown the 
utility of using logarithmic rather than linear binning: much 
of what is interesting in correlations happens at small q, and 
this region is probed better by using logarithmic bins. 



In Figure 2, the results of implementing the MK pre- 
scription are shown, using the GHP max topology and 
the parameterizations (J7J) and (§h . Parameter values used 
for the respective power law and exponential MK points 
were taken from the fit to AA2 of Figure 1. (The effect of 
finite binning was checked by inserting these parameter 
values back into the MK formulae for AA2 and find- 
ing good agreement between UA1 data and MK output.) 
Again, all MK points shown are determined only up to 
an additive constant, so that they may be shifted up and 
down in unison. It is clear, though, that the shape of the 
measured third-order cumulant differs significantly from 
that predicted by the APW formulae and AA2 parame- 
ter values. Similar discrepancies were found when using 
other topologies such as the GHP sum variable || . 
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eral hundred MeV) where Coulomb repulsion is not 
expected to be important. 

• We have checked through additional small-q cuts 
that possible residual track mismatching (multiple 
counting of tracks Q , JL2| ) does not explain the 
discrepancy between predicted and observed AAV 

• Our sample consists of 15% unidentified kaons and 
protons which were treated as pions in our analy- 
sis. This impurity is expected to weaken the corre- 
lations. We have repeated the analysis with a low- 
momentum (|p| < 0.6 GeV/c) sample of particles 
by using the information on dE/dx to get an almost 
pure sample of pions. No reduction of the discrep- 
ancy between the predicted and observed AA 3 was 
observed. 

• We have checked whether the restriction to the azi- 
muthal angle region 45° < \<p\ < 135° distorts AA2 
or A A 3 . In the region q < 1 GeV, no distortion is 
observed. 

Nevertheless, we do not regard our result as a failure 
of the APW formalism per se. One caveat relates to the 
structure of the UAl multiplicity distribution which is 
non-Poissonian (i.e. the cumulants are nonzero at large 
q). The use of an additive constant in the cumulants 
is only an approximate remedy which can be improved 
substantially p3[ . 

Furthermore, we assumed the chaoticity A to be 
momentum-independent and the source currents real- 
valued, as has been done implicitly in all previous ex- 
perimental work known to us. More general versions of 
the APW theory have yet to be tested. Indeed, we hope 
that with our sensitive new testing methods the role of 
such hitherto inaccessible dynamical information can be 
investigated. Beyond such extensions there may lurk a 
non-Gaussian source current. 



FIG. 2. UAl third-order "GHP max" cumulant data, to- 
gether with APW predictions based on Monte Karli integra- 
tion and parameter values taken from AA2- 

This result does not contradict an earlier UAl paper 
|ll| whose higher order moment data was claimed in |5]] 
to be in agreement with the APW formalism. The reason 
is that moments are numerically dominated by a combi- 
natoric background of lower orders, so that comparisons 
based on moments are not very sensitive. Higher-order 
cumulants are a much more sensitive test of the theory, 
as we have just shown. 

We have considered the following sources of errors: 

• No corrections for Coulomb repulsion were in- 
cluded. However, A A3 data rise more strongly 
than theoretical predictions even for large q (sev- 
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